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We measure the field transmission matrix t for microwave radiation propagating through random 
waveguides in the crossover to Anderson localization. From these measurements, we determine the 
dimensionless conductance g and the individual eigenvalues r n of the transmission matrix tt^ whose 
sum equals g. In diffusive samples, the highest eigenvalue, n, is close to unity corresponding to a 
transmission of nearly 100%, while for localized waves, the average of n, is nearly equal to g. We 
find that the spacing between average values of In r n is constant and demonstrate that when surface 
interactions are taken into account it is equal to the inverse of the bare conductance. 
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More than 50 years ago, Anderson [1] showed that be- 
yond a certain threshold in disorder the electron wave- 
function within a material becomes exponentially peaked 
and diffusion ceases. In the intervening years, local- 
ization and its precursors in diffusive samples of en- 
hanced fluctuations and suppressed transmission have 
been shown to affect every aspect of transport [2 -4 . Re- 
lating transmission in disordered systems to random ma- 
trices [5HT2] has provided a powerful approach to calcu- 
lating the scaling and fluctuations of conductance. In 
this approach, the flow of electrical current through a 
disordered conductor is assumed to proceed via a set 
of orthogonal channels on the input and output sides of 
the sample which are coupled via the field transmission 
matrix t. Each channel is a superposition of the N or- 
thogonal transverse momentum channels supported by 
the sample leads. Alternatively, t gives the coupling of 
the field between different points on opposite surfaces 
of the sample. The focus of random matrix theory has 
been the calculation of the dimensionless conductance, g, 
which is the conductance in units of (e 2 //i), and equals 
the sum of the eigenvalues r n of the transmission ma- 
trix tt\ g — X^Li T n- However, direct comparisons with 
measurements of the individual transmission eigenvalues 
have not been made. 

Dorokhov [6 showed that even in conducting sam- 
ples, the current in most channels would be exponentially 
small so that the conductance is dominated by "active" 
or "open" channels with r n > 1/e. The number of such 
channels is close to g, N e ff ~ g [6j [7]. Transmission 
eigenvalues may be expressed in terms of associated local- 
ization lengths, £ n , r n = exp(— L/£ n ) with active chan- 
nels corresponding to the condition L < £ n . Dorokhov 
found that the average spacing between inverse localiza- 
tion lengths of adjacent eigenchannels in a sample made 
up of N parallel chains with weak transverse coupling to 
neighboring chains was constant and equal to the inverse 
of the localization length, l/£. Here, £ = N£ is the local- 
ization length and t is the localization length in a single 
chain and corresponds to the mean free path [5,6]. 



Subsequently, Stone et al. [11 J applied a maximum en- 
tropy hypothesis to random transfer matrices and found 
an expression for r n in terms of eigenparameter A n of 
the transfer matrix or equivalent parameters z/ n , r n = 
1/(1 +A n ) = 1/ cosh 2 (^ n /2). A Coulomb gas model, orig- 
inally introduced by Dyson [13] to visualize the logarith- 
mic repulsion between the eigenvalues of random Hamil- 
tonian, was extended to describe the repulsion between 
"charges" at positions v n . The approximately uniform 
density of charges at v n [lOj [11] suggested a bimodal dis- 
tribution for the density of transmission eigenvalues of r n 
for diffusive waves with peaks near r = and 1 [I2j [14] . 
This is confirmed in microscopic calculations made by 
Nazarov for non-absorbing diffusive samples with perfect 
leads [T5] , 

Recent optical measurements by Vellekoop and Mosk 
have shown that it is possible to enhance transmission by 
manipulating an incident optical wavefront to preferen- 
tially excite high-transmission channels [16 . Total trans- 
mission was enhanced by 44% and focused intensity by 
three orders of magnitude by adjusting the phase of the 
incident optical wavefront reflected from a spatial light 
modulator (SLM) in response to feedback of transmitted 
intensity. Popoff and co-workers [17 used an interfer- 
ence technique to measure a portion of the optical field 
transmission matrix t using an SLM and CCD camera. 
They found that the probability density of the singular 
values of t corresponded to a quarter-circle law [18] [19] . 
The quarter-circle law arises when the elements of the 
t are distributed as independent Gaussian random vari- 
ables and occurs whenever the dimension of measured t is 
smaller than g [20] . Thus these measurements do not re- 
flect correlation in the medium and can not differentiate 
transmission in samples with different values of g. 

In this Letter, we present measurements of the trans- 
mission matrix in the crossover to Anderson localization. 
The transmission matrix is measured for collections of 
randomly positioned dielectric spheres in a copper tube 
for which values of g range from 6.9 to 0.17. The dimen- 
sionality of the measured matrix matched the number 
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of propagation modes of the empty waveguide. For the 
most diffusive sample studied, the average value of t\ is 
0.93, corresponding to an enhancement of transmission 
by 900% over the ensemble average value. The number 
of open eigenchannels with r n /r\ > 1/e with localiza- 
tion lengths greater than the effective sample length is 
equal to g for diffusive waves. For the most strongly 
localized waves studied, the average energy transmitted 
in the highest channel is more than 96% of the conduc- 
tance. We find constant spacing between the averages of 
lnr n in adjacent channels, (lnr n ) — (lnr n +i) = 1/g", and 
identical distributions of In r n for n > 1 in both diffusive 
and localized regimes. We demonstrate that once inter- 
actions at the interface are taken into account in terms of 
effective sample length and localization length, L e ff and 
£, respectively, the constant spacing of values of (lnr n ) 
is equal to the inverse of the bare conductance, giving 
g"=g = £/L e// . 

Because the phenomena we seek to explore depend 
only upon g, the particular sample studied is immate- 
rial. Measurements were carried out on well character- 
ized samples [2F for which we could access values of 
g through the localization transition. The experiment 
setup is sketched in Fig. la. The samples are com- 
posed of alumina spheres with a diameter of 0.95 cm 
and refractive index 3.14 embedded in Styrofoam shells 
and contained within a copper tube of diameter 7.3 cm. 
The field transmission matrix was measured in samples 
with lengths L = 23, 40 and 61 cm over the frequency 
ranges 10-10.24 GHz and 14.70-14.94 GHz with use of 
a vector network analyzer. The wave is localized in the 
lower frequency range and diffusive in higher frequencies. 
The numbers of transverse waveguide modes supported 
in the lower and upper frequency ranges are N ~ 30 and 
N ~ 66, respectively. Measurements are made for two 
linearly polarized components of radiation along the hor- 
izontal and vertical orientations of the source and detec- 
tor antennas between N/2 points on a grid on the input 
and output surfaces of the sample tube in each frequency 
range. New statistically equivalent samples are produced 
by momentarily rotating and vibrating the copper tube 
after the full matrix t is recorded in both frequency ranges 
over a period of two days. The impact of absorption upon 
the statistics of t is removed by compensating the field 
transformed into the time domain for the exponential de- 
cay due to absorption and then transforming back into 
the frequency domain. For the shortest sample studied 
for which g = 6.9, the transmittance is reduced by 25% 
due to absorption. Measurements were made in a total 
of 23 sample configurations for L = 23 cm and 6 config- 
urations for L = 40 and 61 cm. 

Speckle patterns for two different positions a and a' of 
the source antenna in a random realization at a single 
frequency are shown for a diffusive sample with L = 23 
cm, g = 6.9 and for a localized sample with L = 40 
cm, g = 0.37 in Fig. lb. Characteristically, the speckle 




FIG. 1. (Color online) (a) Schematic experiment setup. 
Source and detector wire antennas are mounted on a two- 
dimensional translation stage, so they could move freely on 
a two- dimensional grid. Measurements are made for two or- 
thogonal linear polarizations by rotating the wire antennas, 
(b) Speckle patterns from two different source positions a and 
a at a single frequency in one random realization for diffu- 
sive waves with g — 6.9, (upper row) and localized waves with 
g = 0.37, (lower row). 



patterns for diffusive waves are quite different, while 
the patterns for localized waves are nearly the same. 
The speckle patterns are different for diffusive waves be- 
cause the transmitted field is a random admixture of 
many transmission eigenchannels. In contrast, the ra- 
tio of successive transmission eigenvalues for localized 
waves is large and a single channel typically dominates 
transmission. This is confirmed in the spectra of the 
first few transmission eigenvalues for diffusive and lo- 
calized waves shown in Fig. 2. Also shown in Fig. 2 
is the transmittance, which is the sum of the flux be- 
tween all N 2 pairs of incident and outgoing channels, 
T = ^2ab=i \tab\ 2 - Approximately 7 highly transmitting 
channels with values of r n /r\ > 1/e are seen to contribute 
to the transmittance T with average dimensionless con- 
ductance g = (T) =6.9. In contrast, in the sample with 
g = 0.37, (t2)/(ti) ~ 0.088, and transmission is domi- 
nated by the first eigenchannel for most positions of the 
source. The average value of t\ of 0.93 for diffusive waves 
is nearly a factor of 10 greater than the average value of 
the total transmission for an incident wave from a source 
at point a, (T a ) = (T)/N = g/N ~ 0.1. For localized 
waves, r\ ~ T, so transmission is enhanced over (T a ) by 
a factor nearly equal to N = 30. 

The relative values of conductance for samples of differ- 
ent length and waves at different frequency are given by 
the ratio of the ensemble average of the sum of intensity 
for all pairs of incident and output points and the same 
measurement in a tube emptied of scatterers. The ab- 
solute transmittance is found using a normalizing factor 
obtained by equating (T) = g with the bare conductance 
g = g" = 6.9 in the most diffusive sample. This sample 
is far from the localization threshold and corrections to 
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FIG. 2. (Color online) Spectra of the transmittance T and 
transmission eigenvalues r n for (a) diffusive sample of L=23 
cm with g = 6.9 and (b) localized sample of L=40 cm with 



0.37. The black dashed line gives T = J2a,b=i M 2 
r n and the solid lines are spectra of r n . 



g due to renormalization are expected to be negligible. 
We show below that the scaling g" provides an accurate 
determination of g . 

The statistics of In r n are presented in Fig. 3 for en- 
sembles with different values of g. The probability den- 
sity of individual In r n of the first few eigenchannels and 
their contribution to the overall probability density of 
P(lnr) = Y^n=i ^( mr n) f° r the most diffusive sample of 
g = 6.9 is presented in Fig. 3a. Aside from the fall of the 
probability density P(lnr) near lnr ~ 0, which reflects 
the restriction t\ < 1, P(lnr) is nearly constant with 
ripples spaced by l/g" [22, 23 . We find that the density 
of lnr n for n > 1 is a Gaussian function peaked at (lnr n ) 
with width independent of the index n [24 . The density 
of u = u n = lnr n — (lnr n ) for n between 2 and 10 is 
shown in Fig. 3b and compared with a Gaussian distri- 
bution. The regular spacing of (lnr n ) seen in Fig. 3c 
and the nearly identical distribution of In r n for n > 1 is 
reminiscent of a crystalline lattice arrangement of atoms 
at a finite temperature. 

The slope of the curve in Fig. 3c, is the spacing 
between successive values of (lnr n ), l/g". For diffu- 
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FIG. 3. (Color online) "Crystallization" of transmission 
eigenvalues, (a) Probability density of lnr n (lower curves) and 
the density of lnr (top dash curve), P(\nr) = X)n^( mTn ) 
for diffusive sample with g — 6.9. (b) Probability density of 
u — u n = In T n — (lnr n ) for n between 2 and 10, for the same 
sample as in (a) compared with a Gaussian distribution, (c) 
Variation of (lnr n ) with channel index n for sample lengths 
L = 23 (circle), 40 (square) and 61 (triangle) cm for both dif- 
fusive (green open symbols) and localized (red solid symbols) 
waves fitted, respectively, with black dash lines. 



sive samples, for which there are several channels with 
Tu/ti > 1/e, the number of such channels is equal to g 7 '. 
Dorokhov [6] predicted that g 7 ' is equal to the bare con- 
ductance, g" = g = Nl/L. But the bare conductance g 
should be influenced by wave interactions at the sample 
interface, as is the case for transmission for a single inci- 
dent channel [25H27] . The impact of the interface upon 
transmission can be found by considering the angular de- 
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FIG. 4. (Color online) Illustration of the photon diffusion 
model for the spatial variation of intensity, im- 



pendence of transmission. 

Measurements of total transmission in diffusive sam- 
ples are well described by the expression, T{0) = 
(z p cos# + Zb)/(L + 2zb) = (z p cosO + Zb)/L e ff, where 
is the angle between the normal and the wave as it pen- 
etrates into the sample [27]. This expression is obtained 
from a model in which the incident wave is replaced by 
an isotropic source at a distance of travel from the in- 
terface, z pi at a depth, z p cos9 as illustrated in Fig. 4. 

Diffusion within the sample gives a constant gradient 
of intensity to the right and left of the effective random- 
ization depth z p cos 0, which extrapolates to zero at a 
length z\y beyond the sample boundary. The gradients 
of intensity at the input and output surfaces give the 
reflection and transmission coefficients, which together 
with the condition, T + R = 1, gives the expression for 
T(ff) above. The linear falloff of intensity near the sur- 
face boundaries should hold even for localized waves since 
the diffusion coefficient varies with depth into the sam- 
ple but is hardly renormalized near the sample bound- 
aries [28 . Since z v and z^ are proportional to the mean 
free path, averaging over all incident angles gives, 
g = r]N£/ (L + 2z h ) = riN£/L eff = £/L eff . Here rj is 
independent of L and includes the effects of reduced flux 
into the sample due to external reflection of the incident 
wave at the interface, and enhanced internal reflection. 
These effects tend to cancel in transmission so that rj ~ 1. 
The values of z^ of 13 and 6 cm for the diffusive and lo- 
calized samples, respectively, are obtained by fitting the 
diffusion model to the measured time of flight distribu- 
tion [29, 30 . Hypothesizing that g" — g , the product 
g['L e ff should give a constant length, which we can iden- 
tify as £. We find this product to be essentially constant 
for measurements carried out in the three sample lengths 
in each of the two frequency ranges as shown in Fig. 5. 

The nearly uniform density of P(lnr), seen in Fig. 3a, 
which is associated with the equal spacing between suc- 
cessive values of (lnr n ) is consistent with calculations of 
Dorokhov [SJE]. These results correspond to a nearly uni- 
form probability density P(lnr) = g for diffusive waves, 
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FIG. 5. Identification of g" with the bare conductance g . 
The constant products of g"L e ff for three different lengths 
for both diffusive and localized samples give the localization 
length £ in the two frequency ranges. 



apart from a small modulation due to the peaking of this 
density at lnr = (lnr n ). This leads to the probabil- 
ity density, P(r) = P(lnr)^^ = g/r. This distribu- 
tion has a single peak at low values of r in contrast to 
later prediction of a bimodal distribution for r for dif- 
fusive waves which has a second peak near unity. This 
difference may be a consequence of differences in the cir- 
cumstances of our measurements as compared to the as- 
sumptions made in many analytical calculations. In our 
measurements, the wave is a vector, the sample is not in 
the extreme diffusive limit of g ^> 1, t couples measure- 
ments on a grid of points, and surface interactions are 
included. Another difference between measurements and 
calculations is that flux is not conserved in our sample 
so that the measured transmission matrix is not Hermi- 
tian. However, the spacing of transmission eigenvalues 
does not appear to be affected appreciably by moderate 
absorption. In the most diffusive sample studied with 
g = 6.9, the spacing between adjacent values of (lnr n ) is 
uniform and remains so when the field spectra are com- 
pensated for absorption. Although transmission is sup- 
pressed by 25% by absorption, the average spacing is 
only reduced by ~ 2%. The insensitivity of the statis- 
tics of transmission eigenvalues to absorption parallels 
the insensitivity to absorption of the variance of total 
transmission normalized by its ensemble average, which 
is suppressed by ~ 1.6% [2T] , 

In conclusion, we found the transmission eigenvalues r n 
and the optical transmittance T for microwave radiation 
propagating through random waveguides for diffusive and 
localized waves. We show that the scaling of (r n ) and g is 
determined by the bare conductance, g = £/L e ff, which 
is strongly affected by surface interactions. These mea- 
surements make it possible to directly test random ma- 
trix theories based on different assumptions and provide 
a systematic basis for controlling the transmitted field 
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for applications in communications, imaging and lasing 
in complex systems. 
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